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Abstract 

A representation  of  controllable  linear  systems  is  introduced,  which  permits  assigning  poles  or  charac- 
teristic parameters  to  a state  feedback  system  by  a matrix  multiplication.  This  is  used  as  a link  between 
state  space  and  classical  parameter  plane  methods.  The  system  representation  maps  a point  in  a n*p  dimen- 
sional parameter  space  9 of  characteristic  parameters  into  the  nxp  dimensional  parameter  space  \ of  state 
feedback  gains,  where  p is  the  number  of  actuators.  For  p*l  the  coordinates  of  the  9 space  are  the  coef- 
ficients of  the  closed  loop  characteristic  polynomial,  for  p>l  they  are  coefficients  in  a characteristic 
polynomial  matrix  and  its  determinant  is  the  characteristic  polynomial.  3y  this  computationally  simple 
napping  procedure  it  becomes  feasible  to  map  not  only  a fixed  set  of  eigenvalues  but  also  regions  in  the  s 
or  z plane,  in  which  the  eigenvalues  shall  be  located.  This  relaxation  of  the  dynamic  specifications 
permits  satisfying  other  typical  design  specifications  like  robustness  with  respect  to  sensor  and  actuator 
failures,  .arge  parameter  variations,  tinite  wordlength  implementation,  and  actuator  constraints.  All 
tradeoffs  between  such  requirements  can  be  made  in  the  X space.  Three  examples  illustrate  the  variety  of 
problems  which  can  be  tackled  with  this  new  tool. 


INTRODUCTION 


Control  system  specifications  ire  usually  not  given 
in  terms  of  a quadratic  cost  function  or  i set  of 
eigenvalues.  These  are  mainly  used  as  free  para- 
meters in  trial  and  error  design  procedures  aimed 
at  good  tradeoffs  between  the  dynamics  of  the 
system  and  other  design  aspects,  for  example 
actuator  limitations  and  robustness  with  respect 
to  sensor  or  actuator  failures  or  other  large 
parameter  variations.  Three  questions  in  this 
context : 

1.  Quadruplex  techniques  (for  example  in  aircraft 
control  systems'  are  an  expensive  solution  to 
the  reliability  problem.  Is  it  not  sufficient 
to  guarantee  chjc  ill  unstable  ar.d  insuffi- 
ciently damped  modes  remain  observable  and 
controllable  under  all  considered  combinations 
of  sensor  and  actuator  failures? 

2.  Which  system  cnanges  are  so  essential  that 
they  require  adaptation  of  the  control  system 
by  identification,  failure  detection  etc.? 
Which  range  of  such  changes  can  be  covered 
satisfactorily  by  fixed  gain  feedback  or  a few 
secs  of  gams  and  a simple  switching  criterion. 

3.  Is  a given  set  of  mixed  specifications,  e.g. 

cn  bandwidth  and  damping,  ictuator  constraints, 
robustness  requirements , etc.,  compatible 
within  an  assumed  control  system  structure,  or 
which  of  the  specifications  are  conflicting, 
how  far  do  we  nave  to  relax  them? 

This  paper  does  not  give  a complete  answer  to  these 
questions,  however  l method  is  proposed  and  some 
C jo Is  arc  provided  t j attack  such  questions  under 
the  following  assumptions. 

1.  Only  linear  plants 

i’ii-la.  &■  l*j.  • •*„]  ’ - im  iuj- • 1 (i) 

or  x(kfL)  *A”(k)  ♦Bui*)  are  considered.  It 
is  assumed  chat  oq . (1)  Is  -r it ten  in  "sensor 
coordinates , " i.e.  ail  measured  vjrl.toles  ire 
state  variables  x^.  Sovora*  rairs  lAi»3i), 
(A^.S^),  etc.  may  be  given,  j.g.  for  different 
operating  points  of  an  underlying  nonlinear 
system.  It  is  assured  that  ill  pairs  (A^,3X) 
are  controllable  and  have  tn*  »ane  control- 
lability indices. 


This  research  wa*  supported  by  the  Deutsche  Forscnungs -und 
the  U.  S.  Air  Force  under  Crant  AfOSR  78-36)3. 


2 # 

The  price  of  an  actuator 
with 

is  assumed 

Co  increase 

IT  - max|ui(t)| 

and/or 

(2) 

u ' * max!  (t ) | 

(3) 

wnere  the  worst  initial  state  within  given  limi- 
tations is  considered,  u^  and/or  u?  should  be 
kept  "small." 


3.  A stace  feedback  structure 

u ■ -K  1 £ (4) 

is  assumed.  The  details  and  examples  are  worked 
out  for  single  input  plants  with 

u * -k'x  - -[k1  k,  ...  kjx.  (5) 

For  multi-input  plants  the  basic  result  is 
stated  in  the  Appendix.  The  nxp  elements  of  k' 
are  the  free  parameters  of  the  proposed  method. 
They  are  coordinates  of  a parameter  space 
called  "state  feedback  gain  space"  or  "Y-space." 

U.  It  is  assumed  that  sensor  and  actuator  failures 
occur  in  the  form  that  the  output  of  a failed 
element  is  not  correlated  to  its  tnnut.  Then 
the  output  is  an  external  disturbance.  Rejec- 
tion of  external  disturbances  is  net  considered 
in  cnis  paper.  For  the  closed  loop  system 

x » fA-B  K'  )x  (6) 

in  sensor  coordinates  a failure  of  the  ich 
sensor  (actuator)  is  equivalent  to  a structural 
change  by  which  the  ith  column  (row)  of  K* 
becomes  zero.  Also  in  the  case  that  a state 
variable  is  not  measured,  the  corresponding 
column  of  K1  Is  zero.  It  is  part  of  the  design 
to  decide, which  state  variables  »rc  measured 
and  for  which  of  them  redundancy  must  be  pro- 
vided for  high  svsten  reliability.  It  is  a 
goal  to  avoid  failure  detection  and  multiplexed 
components  whenever  possible. 

3.  It  is  assumed  thit  the  nom'.tial  dynamic  heuavior 
of  the  control  system  can  be  specified  by  a 
region  In  the  eigenvalue  plane  -e.g.  for  a con- 
tinuous time  system  in  the  s plane  the  region 
to  the  left  of  the  bour.dar**  marked  with  3*i  in 
Fig.  I - where  all  eigenvalues  must  oe  located 
for  the  nominal  structure  and  oararoter  values. 
For  failure  situations  a relaxed  "emergency 
specification,"  e.g.  the  boundary  P*0.5  or  th« 
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Fig.  i.  A family  of  hyperbolic  boundaries  in  s 
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stability  iitr.it  -=*0  in  Fig.  I cay  be  given. 

A specification  is  robust  under  a failure  if 
no  eigenvalue  crosses  the  corresponding 
boundary  due  to  the  failure. 

The  proposed  method  is:  Design  inY-space.  \ 
region  Yl  in  the  Y-space  is  determined,  such  that 
all  eigenvalues  meet  the  speci f ications  iff  K'^Yp 
This  ray  be  the  intersection  of  several  such 
regions  for  different  parameter  values.  Subsets 
of  X\  with  certain  robustness  properties  can  be 
found  and  tradeoffs  with  actuator  constraints, 
bandwidth  requirements,  etc.  can  be  made  in  the  Y* 
space.  The  method  also  shows,  whether  a solution 
exists  under  the  given  assumptions  and  if  not, 
which  alternatives  exist  for  relaxation  of  speci- 
fications such  that  a solution  will  exist. 

Parameter  space  methods  have  a long  tradition, 
mainly  in  Russia  and  Yugoslavia.  Siljak  (Ij  gives 
a historical  review  of  the  work  by  Vishnegradsky , 
Neimark,  Mitrovic,  and  others.  Siljak  generalized 
these  parameter  mapping  methods  significantly.  A 
typical  procedure  for  a continuous  time  system  is 
to  assume  a controller  structure  with  two  free 


parameters  o and  3.  Determine  the  closed-loop 
characteristic  polynomial 


?(s)  - 

n 

iS0? 

L(z,i)sl  - 0. 

(7) 

Subs 

Citute  s = r+j 

u and 

separate  eq.  (7) 

into 

its 

rea  1 

and 

imaginar 

7 ?ar 

c s : Rc  (t  ,x',^,3  )* 

0, 

Ia(- 

, r , u , 

,3)*0.  Assume 

Ci-.cse  nonlinear 

cauat ions 

have 

a solution 

o ■ 

o(C,x), 

3 

- 3 (r,*). 

(3) 

Equal 

:ion 

(8)  alio- 

VS  CO 

map  : pairs  on 

the 

boundary 

into  the 

i-3- 

olane.  The  image 

bound 

laries 

divide  the  ^-3-?lor.e  into  regions  characterized  by 
the  nunbor  of  eigenvalues  inside  and  outside  the 
s*plane  region. 

In  the  present  paper  the  control  system  structure 
is  restricted  to  state  feedback.  This  permits 
simplifying  the  determination  of  cq.  (S^  by  pole 
placement  methods.  ConsiJer  for  example  1 second 
order  single  -input  system  with  k^  = v,  k*>»3  in  eq. 
(5).  In  classical  parameter  ol.me  methods  Pfs)  * 
dot  (s  k')  »o  fs  ,5)*T>- (},3  ) is  determined 

and  with  s*?-*-ji  solved  tor  > and  3.  In  the  method 


terms  of  * anJ  u by 

P(s  ) » (S-7+juj)  (s-7-j  L’)  * S“-2?S*7‘4«i T 
-P0(T,.JJ)  +p1(7)s+S2  *0. 

Then  by  pole  placement 


(9) 


kl  ’ ^p0'pl>  ” 


k2  ’ 3(po*Pl)  * 3(?,ji). 


(10) 


Thus  Che  mapping  equation  (8)  is  obtained  in  a 
different  way.  More  generally  for  an  nth  order 
single  input  system  in  both  approaches  an  n dimen* 


sional  parameter  spaced  with  coordinates  p^  Is 


introduced  as  an  intermediate  step  between  the  set 
of  eigenvalues  A*i\j  ...\n]  and  the  ?C-space.  The 
relation  between  A and  X can  be  expressed  in  boch 
directions  : 


a) 


b) 


From  X to  S by  the  characteristic  equation 
P(V  )“det  0\  I-A+bk ' ) , from  9 to  A by  numerical 
factorization  of  pp). 

From  A to  S by  multiplication  of  elementary 
factors  P(\  )" (V  i>(\ -A 2 ) • • • Q *\ n)  > from  $ t0 
X by  pole  placement. 


In  the  nexc  seccion  pole  placement  is  formulated 
in  a form  which  makes  direction  b)  far  more  attrac- 
tive than  direction  a).  In  Section  3 the  use  of 
such  boundaries  for  the  design  of  robust  control 
systems  in  . Y-space  is  discussed.  Section  4 shows 
Che  application  to  three  examples. 


3QUMDARY  MAPPING 


Po  lo-p  laceirent  for  sine  !e-inout  systems.  The  pole 
placement  theorem  is  used  in  the  form  originally 
published  in  German  in  (2),  available  in  English  in 
[3]  : Given  an  nth  order  monic  polynomial  PP. ),  an 
nxn  matrix  A and  an  n/1  vector  b such  thac  decR^O, 

R ,(b,Ab An'',b]  , the  equation  ?p  ) * 

detp.^-Amb  k1 ) has  a unique  solution  and  this  solu- 
tion is 

k1  - e’P(A) 

where  e*  is  the  last  row  of  R*l. 


(11) 


• +)  ‘ 


(12) 

(13) 


With  Pp.)  - (X-Xl)(V-\2)  •••  (X-'„) 

■po+PlX  +-"+?n-lXn'1' 
eq.  (11)  may  be  written  as 

k ' - c ' (A-\  tI)  (A-*  ,T)  ...  (A-'  nI) 

■«'  (P0I  + plA  + • • • +pa-idn  l+An)- 
For  mapping  from  9 space  to  V space  it  is  more  con- 
ventenc  to  rewrite  eq . (15)  as 


(14) 

(15) 


■ o’E 


ri’  1 

d A | 


,[po  p 


Kn-1 


(16) 


• 

L-  -J 


E_  is  an  (n+l)vn  matrix.  If  the  inverse  of  eq . (16) 
is  needed,  it  is  convenient  to  express  the  last 
row  of  £ by  the  Cay loy-Hami Iton  theorem  in  terms 
of  the  previous  rows.  This  however  requires  the 
evaluation  of  the  characteristic  polynomial 


act  0 I -A)  •*  a . * + • • • + 1 


tible  form  ot  eq.  (16)  is 


n-l 


i n-ia 


The  inver- 


1 ^ Po"au 
O 0 


p,-a 


l “l 


pn-l*; 


-1 


n-f - 


e 'A 


! 'A 


n-l 


(17) 


The  columns  of  V can  be  evaluated  recursively  by 
Leverncr's  algorithm,  which  also  gives  the  [2). 


\ . 


1 


! ! 


i ' 


\) 

! 


» 

1 

: 

I 


The  (ip  (16)  is  pose  convenient  foe  trnl  and 
error  design  procedures,  graphical  displays  of 
cross-sections  of  the  v space,  etc.  The  plant 
description  in  the  form  of  the  matrix  E is  evalu- 
ated only  once  for  a given  pair  (A,b).  The  map- 
ping of  a trial  design  point  in  ->^space  then 
requires  n~  multiplications  and  n-  additions.  This 
compares  favorably  with  mapping  a trial  Jesign 
point  from  the  parameter  space  of  quadratic 
criteria  via  the  Riccati  equation  into  \ space. 

The  generalization  of  eq.  (16)  to  multiinput 
systems  is  given  in  the  Appendix. 


Sensitivities . The  influence  of  a coefficient  p^ 
of  the  characteristic  polynomial  on  k ' , given  the 
other  p.,  follows  from  eq.  (16)  as 

J d-  > .'A1 

377  - - ■ (l8) 

The  influence  of  an  eigenvalue  on  k ' , given  the 
other  'j,  follows  from  eq.  (14)  as 

(A-\  jl)  . . . (A j.jl)  <A-\  1+1I)  . . . (A-\„I)  , 


'(19) 


For  complex  conjugate  eigenvalues  quadratic  factors 
in  P(s)  are  more  convenient.  Let  P(\  ) s 
(a+b\-»\2  )q(\  ) f then 


e’  (al  +bA  + A2)-Q(A) 

. d k * , 

e’Q(A),  -Tr-  = e'AQ(A). 


(20) 

(21) 


Rerrions  in  2 olanc . Boundaries,  symmetric  with 
respect  to  the  real  axis,  in  the  ' plane  will  be 
specified,  which  describe  the  desired  eigenvalue 
locations.  Two  cases  will  be  discussed:  A real 
rooc  crossing  a boundary  or  a complex  conjugate 
pair  crossing  a boundary.  The  third  possible  case 
of  roots  leaving  the  region  through  infinity  can 
be  avoided  by  closing  the  contour  by  an  arc  of  a 
circle  with  large  radius.  Typically  the  region  is 
a connected  set  and  the  boundary  has  two  inter- 
sections with  the  real  axis.  In  this  case  there 
are  two  real  root  boundaries  and  one  complex  root 
boundary  in  X space.  However  Dther  boundaries  are 
possible,  e.g.  separate  boundaries  for  slow  and 
fast  modes,  etc.  For  each  intersection  of  a 
boundary  in  the  \ plane  with  the  real  axis  at  \ =*7 , 
a real  root  boundary  in  X space  is  obtained  from 


?0.)  - (X-r)R(X),  R0.)sr  +:,• 
o i 


+rn  .X"-2*"-1. 

(22) 


By  eqs . (16)  and  (22) 

WPo  •••  Vl>  “ki(--ro  •••  V2>  i*1’2’-'1  (23) 

where  the  k*  depend  linearly  on  ? . Thus  7 can  be 
eliminated  by  one  of  the  k^'s  to  give  the  linear 
boundary 

k3  •••  rn-l)  J’1-2 n-  j*1-  <24> 

This  is  a straight  lino  in  tho  k : - k . plane. 

Another  part  of  the  boundary  is  obtiined  if  a 
complex,  pair  of  eigenvalues  crosses  the  boundary 
at  \ ■ r-t-jjj.  Then 


?0  ) - (\‘-2r\  +t- 
Q0 ) -q*  -nq, 


x*)Q(\) 

Cq.  _\n"3  +'  n~2 


(25) 


''o  'U'  Mn-r 

ind  the  type  of  boundary  in  X space  depends  on  the 
form  of  the  boundary  r-j.(r)  in  the  \ plane.  For 
: "const.,  i.e.  a parallel  to  the  imaginary  axis, 
Che  k^  arc  linear  in  e* . Thus  for  given  the 
image  in  the  X space  is  a straight  line. 


Most  cormonly  used  boundaries  are  conic  sections 
symmetric  to  the  reil  axis,  i.e. 


+c  - . 

Special  cases  arc  * 


cn < 0 ellipse,  of  particular  interest  are  circles 

c 7 * - 1 , e.g.  constant  natural  frequency  curves 
in  s plane,  stability  limit  and  other 
boundaries  in  z plane. 

Co-O  parabola,  or  if  also  c^O,  c >0  straight  line 
parallel  to  the  real  axis. 
c1>0  hvpcrbola,  in  particular  2 straight  lines  for 
i-"C2 (J-^o)4* » » e.g.  constant  damping 

lines  in  s-plane. 

Figure  l shows  the  family  of  hyperbolas 
2 2 2 2 

u)  » -P  +7T/P  . (27) 

For  c-0  this  goes  to  the  imaginary  axis,  for  o«l 
the  asymptotes  are  the  l/v‘2  damping  lines.  For  a 
different  scaling  z may  be  replaced  by  7/d. 
Substituting  eq.  (26)  into  (25) 

P(X)  - [>.2-2<7\  + (l+c2)a2  + cf  +co)QG),  (28) 

This  shows  that  the  and  depend  linearly  on  e 
for  c2=-l,  i.e.  for  circular  boundaries.  Thus  also 
in  this  case  the  boundaries  in  X space  are  straight 
lines  if  the  eigenvalues  in  Q 0. ) are  fixed  and  a 
complex  pair  of  poles  moves  along  the  circle. 

For  Cjj1-!  the  functions 

?i  ■ pi(c-q0  •••  Vj>  (29) 

are  quadratic  in  ? and  by  eq.  (16)  the  same  is 
true  for 

k.  - k.(c.q0  ...  qn.3).  (30) 

One  of  these  equations  can  be  solved  for  7“?(k.), 
where  only  real  roots  7-v^5  are  of  interest,  if 
only  the  left  half  \ plane  branch  of  the  conic 
section  is  needed,  then  7=^-6  is  selected  and  sub- 
stituted into  the  other  k^  equations  to  give 

-kj(ki,qo  ...  qn.3),  i*j.  (31) 

Note  that  this  is  not  the  curve  in  the  k. -k.- 
cross- section  of  the  X space,  which  would  hi 
obtained  for  k^=const.f  mfi,j;  eq.  (31)  gives  the 
curve  kj(k^)  tor  constant  q . ..qn-V  For  the 
numerical  determination  of  boundaries  In  cross- 
sections  of  the  'V  space  the  implicit  form  (30)  is 
more  useful.  It  gives  7-values  as  a parameter 
along  the  boundary.  Constant  damping  spirals  in 
the  z plane  are  not  conic  sections.  Usually  they 
arc  supplemented  by  a condition  \ z\  <a,  a<l.  The 
resulting  regions  can  be  reasonably  well  approxi- 
mated by  a family  of  nonintersecting  circles 
2 2 2 

(v-Vq)  +v  * r , z ■ v + jw 

V (v  -1)  - 0. 99r (r-l)  , v < 0.5.  (32) 

o < J o 

It  is  shown  in  Fig.  2.  For  r«l  it  is  the  unit 
circle, with  decreasing  r the  center  v of  the 
circles  moves  to  the  right  until  it  reaches  0.45 
for  r=0.5,  it  then  goes  back  to  zero,  where  v mr"0 
is  the  deadbeat  solution. 

Regions  in  » space.  Equations  (22)  and  (25)  show 
that  the  mapped  boundaries  in  V space  represent  the 
conditions  under  which  the  number  of  eigenvalues 
inside  and  outside  a \ -region  can  change.  These  V- 
boundarics  partition  the  v space  into  regions;  each 
of  them  corresponds  to  a fixed  number  of  eigen- 
values inside  the  \ region,  and  it  must  be  decided, 
for  which  V region  all  eigenvalues  are  inside  the  \ 
region.  For  closed  contours  in  Che  '-plane  the  v 
region  is  bounced,  since  by  eq.  (16)  no  kj,  can  go 
to  infinity.  If  there  are  several  bounded  regions, 
a simple  test  is  to  check  the  eigenvalues  for  an 
arbitrary  K'  in  the  considered  V region.  An  alter- 
native are  Sitjak's  "shading  rules"  for  the 
boundari*-  i ll . 

Consider  a second  >rder  system  and  a circular 
region.  Boundaries  in  the  kj-kn-piane  are  three 


i 
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i. 


(26) 


Fig.  2.  A family  of  circular  boundaries  in  z plane. 

straight  lines  obtained  for  the  two  real  root  cases 
' and  \»r->  and  one  complex  root  case.  They  par- 
tition the  v plane  into  seven  regions  with  the  pro- 
perties: no  pole  outside  the  circle,  one  left,  two 
left, one  right,  two  rignt,  one  left  and  one  right, 
complex  outside.  The  only  bounded  region  is  the 
triangle,  thus  for  k^,’^  in  the  triangle,  \ and 
are  inside  the  circle.  At  the  vertices  of  the  tri- 
angle both  poles  cross  boundaries  simultaneously. 
This  is  Che  case  for  L)  a double  pole  at  7,,  2)  a 
double  pole  at  ^2 , and  3)  one  pole  at  and  one  at 
~2  . Thus  the  total  rapping  procedure  consists  of 
jusc  three  pole  placements,  i.e.  twelve  multiplica- 
tions and  twelve  additions.  This  makes  it  easy  to 
map  the  family  of  circles  of  Fig.  2.  If  the  circle 
is  deformed  to  a different  closed  contour  with  the 
sane  real  axis  intersections  at  7 ^ and  7^,  then  the 
three  vertices  and  two  edges  of  the  triangle  remain 
unchanged,  the  third  edge,  i.e.  the  complex  root 
boundary,  is  replaced  by  a curve.  For  a third 
order  system  and  a circular  ' -region  the  two  real 
root  boundaries  are  planes.  The  complex  root 
boundary  for  a fixed  real  pole  in  a straight  line. 
3y  moving  the  real  pole  the  straight  line  moves  and 
forms  the  third  surface.  The  4 vertices  of  the 
region  arc  again  obtained  by  pole  placement  of  the 
four  characteristic  polynomials  with  ceros  in  the 
sec  For  the  correspond ing  region  in  •? 

space,  Fam  and  Med  itch  (a)  have  shown  chat  the 
convex  hull  of  the  region  is  the  tetrahedron  with 
the  above  mentioned  vertices.  3y  the  linearity  of 
the  turning  equation  k’^n'F  this  property  does  hold 
in  the  v space  also.  Similarly  f >r  arbitrary  n 
from  [4]  follows:  The  convex  hull  of  the  v region, 
f^r  which  all  eigenvalues  are  located  inside  a 
circle  with  center  on  the  real  axis  and  inter- 
secting the  real  axis  <it  is  a polyhedron 

with  n-rl  vertices.  They  can  5e  obtained  bv  pole 
placement  tor  the  r.-rl  cnaracter iscic  polynomials 
with  ceros  in  the  set  1 7 p 7 ^ ' • The  two  real  root 
hypernlanes  are  two  of  the  surfaces  of  the  poly- 
hedron. 

3.  of,si c.i  in  v s?.\r.r. 

Robuscn^si  with  t>  -.erv.  ’r  : u ir-a  . A 

specification  is  "F,  -r>busc"  if  it  remains  satis- 
fied after  a failure  >f  sensor  i,  it  is  "F^:- 
robust”  if  the  same  holds  after  failures  of  both 
sensors  1 and  j.  Fig.  3 shows  an  example  of 


Fig.  3.  Illustration  of  failure  robustness  and 
emergency  boundaries. 

boundaries  in  the  kpk0-subspace . It  is  assumed, 
that  all  other  feedback  gams  are  fixed.  The  pro- 
jection of  point  1 on  the  k,  axis  is  outside  the 
energency  region,  i.e.  the  emergency  specification 
is  not  F2*robust,  it  is  however  F,  robust.  Points  2 
and  3 are  F^  and  F2  robust.  No  point  is  Fp  robust 
since  the  origin  k,*k^*0  lies  outside  the  emergency 
region.  Point  3 also  meets  the  nominal  specifica- 
tion and  is  a good  candidate  for  a robust  control 
system.  Since  the  nominal  boundary  intersects  the 
ko  axis,  an  alternative  to  the  robust  solution  3 is 
to  eliminate  the  x 1 sensor  and  to  multiplex  the  x, 
sensor.  This  voulu  maintain  the  nominal  specifica- 
tions under  a failure  ofone  of  the  x*»  sensors. 
However  it  requires  failure  detection  with  at  least 
three  x~  sensors  and  very  likely  is  the  more  expen- 
sive solution. 

Robustness  with  rcsnect  to  actuator  failures. 

Assume  kj  and  k n in  Fig.  3 are  elements  of  different 
rows  of  the  feedback  matrix  K ' in  eq.  (A.S).  Then 
the  same  arguments  as  above  apply  for  the  robust- 
ness of  specifications  with  respect  to  actuator 
failures . 

Robustness  with  respect  : > iarge  parameter  v.u’.j- 
t ions . Assume  that  for  different  operating  condi- 
tions different  pairs  (Ai^)  etc.  are 

given.  One  boundary  in  the  eigenvalue  plane  now 
maps  into  different  boundaries  in  \ space,  and  it 
must  be  tested,  whether  there  exists  an  intersec- 
tion of  the  admissible  regions.  ff  it  does  not 
exist,  then  at  least  1 gun  scheduling  system  can 
be  designed,  in  which  each  gain  covers  as  many  oper- 
ating conditions  as  possible. 

R •‘hur.  vifh  t • f:n:t  • >r  ! I or-- 1 •; . The 

feedback  control  law  may  be  implemented  approxi- 
mately in  a short  wordlength  micr oprccc s sor  as 

u-hlu  ■ (K'+^K*  ) (X'+ix)  a*  ,'K'x  + K*  (33) 

For  small  x Che  dominant  term  m in  is  K'lx,  I.e. 
the  gains  should  be  not  Coo  high.  For  large  the 
dominant  term  is  \K*v . The  maximally  IV.'  r^oust 
solution  is  the  center  of  the  largest  hyoercube 
with  edges  parallel  to  the  axes  in  the  admissible 
Y-rcgion.  Fig.  * illustrates  „\K'  robustness. 

Ac  t'l.i ter  constraint-,.  Constraints  on  u*"\ax|u(t>! 
and  u»max|  On  t)i  can  be  treated  in  Y space.  For  the 
regulator  problem 

I u (C ) I • !k\(e)|  i k.  X (C ) (34) 

•Jith  equalicy  for  Che  -.or-iC  ;ase  of  x(t)  ie.g., 

<_*ck  for  sone  c*0).  Assuming  chjt  all  state  vari- 
ables have  beer,  nor-jiiced  Co  their  naxiiauo  value, 


I 


Che  norm 


with  damping  l/\  2 or  more  is  required: 
Boundaries  A and  3 in  the  s-plane  of  Fig.  5. 
This  specification  shall  be  Fo  , F3 , and  F03 
robust. 

iii)  A tradeoff  with  the  magnitude  of  u*max|u|  and 
the  maximum  bandwidth  is  to  be  ;»uie. 

The  observability  analysis  shows , that  j cannot  be 
observed  by  1 or  i,  thus  the  j sensor  is  essential 
i.e.  the  reliability  with  respect  to  a failure  of 
the  y sensor  can  be  increased  only  by  redundant  y 
sensors.  The  gain  will  be  determined  first.  r-> 
robustness  of  stability  requires  chat  P(s)  ■ 
s +(1-H: )s2+(1-k: )s+k^ , c>0,  is  Hurvitz,  i.e. 

0<k|<  (i+c)-.  The  worst  case  is  c— 0.  For  maxi- 
mum bandwidth  choose  ki"i. 


i.e.,  Che  distance  from  the  origin  in  V space  can 
be  used  as  a measure  for  u~max i u ( t ) | . Similarly 
|u(t)|  ■ik'x(t)|  * Ik' (i*b  k' )x(c)| 
and  k (A-b  k ' ) can  be  used  as  a measure. 


F.XAltPT.ES 


The  following  three  examples  of  second,  third,  and 
fourth  order  show  various  typical  design  aspects 
and  solutions  in  v space  using  the  tools  intro- 
duced in  this  paper.  All  calculations  were  done  on 
a programmable  pocket  calculator. 

Second  order  discrete  svsten. 


x_(k+l ) * A x (k)+  ou  (k) , A 


Find  u*-lk^  k0 1 x such  that 
1)  stability  is  F^-robust  and 
?2 -robust , 

ii)  the  system  remains  stable 
with  the  maximum  _k. 


possible,  also 


i.e.  kj-pQ-l,  ko*P| -cp0 - l+c*- , k^*p0-l-c.  By  eq. 

(25)  P(s  Ws2-2  )+(c*--K2-27q)s 

+ ( q-2r  )s^+s^  . P0*q  (•T'4,i-)*1  guarantees  q>0,  i.e. 

Che  real  eigenvalue  is  stable.  For  c *2  and  boundary 
A,  i.e.  1" , P -j " 2 qr “ , p»*2r(r-q),  P2*q-2~  and  by 
eq.  (41)  k3*2qc““l , k;»2::  (r+2  )-2q  (l+r  )+3 , k3-q-2r-3. 
3y  kj*p0*l,  q=l/2'-,  Che  product  of  eigenvalues  is 
fixed  and  k.,*2?  (c+2)- (l+r  )/?-*3  t ld«l/2r' -2c-i  . 

This  is  curve  A in  Fig.  5.  A Bode  ploc  shows 


The  vertices  of  the  stability  triangle  in  the 
k^-ItT-Flane  are  determined  by  2 pole  placec.encs 

1.  P(a)  - (:-*l)2  . k ' “[21  6] 

2.  ?(s)  * (z+l)(z-l)  * z:-l,  k ’ = [ - 1 -14]  (38) 

3.  ?(s)  - (1-1)4  .-2.2-+1.,  k ' =1-3  - 10)  . 

In  Fig.  - it  is  seen  that  the  requirements  for  F| 
and  F-,  robustness  are  not  compatible,  the  F ^ robust 
region  is  chosen.  ii)  requires  to  place  the 
largest  square  with  sides  parallel  to  the  axes  into 
the  r ^ robust  region.  It  has  the  center 
k ' *[  -o . - 5-0-0  -10.7272721  and  permits  Ik*  1 . 4545^5 . 

This  k'  places  the  eigenvalues  at  =0.132, 

Z'y  *0.636. 


s P'cne 


z2 , Pj  ana 
Poousf  *^eq.cn 


r2  PcCuSt 


Max  Bandwidth 

Fig.  5.  Tradeoffs  between  tachometer  (lu)  and 

amperemeters  (k3)  feedback  for  a dc  motor. 

that  the  maximum  banJwidth  is  obtained  if  real  and 
complex  poles  have  the  some  distance  from  the 
origin,  i.e.  q*l,  :*-l/v2,  see  Fig.  5.  The  minimum 
u is  obtained  at  the  point  of  curve  A,  which  is 
closest  to  the  origin. 

For  boundary  B P(s  )*  (s -?)*■  (s+q  )*qc--t-7  (?-2q  )s  + 

(q-2r )s-+s^  , and  as  before  ko®"'  (r+4) -2  (l*a )/?*+3 , 
ky3 l/cfc -2r-3 . Points  in  the  shaded  region 
satisfy  specifications  ii).  The  point  of  maximum 
bandwidth  on  the  A curve  is  nothin  the  robust 
region,  therefore  damping  > 1/s  2 is  necessary.  If 
q*l  is  kept  constant  ind  the  complex  poles  move 
along  the  unit  circle,  then  k moves  along  a 
straight  lino  in  the  k2 -kj-p lane . Two  points  of 
this  line  are  k-*“k3*0  for  a triple  pole  at  s*-l  and 
the  point  ~ ■-l/\T  on  curve  A.  Thus  a good  first 
choice  would  be  the  point  k*>“O.J,  kj^-O.S  indicated 
by  the  triangle  in  Fig.  5. 


’ODjs: 


,k  - PoCust 


Fig.  4.  Second  order  system:  A circle  maps  into  a 
triangle.  Finite  word  length  robustness. 

DC  servo  motor.  State  variables  in  sensor  coordi- 
nates are  x*l :»  i i 1 * with  >*  shaft  angle,  •'  ■ 
angular  velocity,  i "armature  current,  input:  volt- 
voltage  u.  Assumptions  : a)  load  torque  >!l*c;c, 
b)  all  state  variables  normalized  to  their  maximal 
values,  c)  simple  numerical  values 


d)  scace  feedback  u*  k,  (r-xj^-k^x-t-i^xj , (40 

which  gives  a zero  stationary  error  l imi  r-\i (t ) ) 
for  a step  reference  input  r,  prwiA-d  the  system 
is  stable.  Find  k’"ik^  kj  kjl  such  that 

i)  stability  la  r-13  robust  for  all  1 'ads  c. 

ii)  For  a load  c*2  i complex  pair  of  eigenvalues 
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C r i n.» . Consider  n crane  with  the  physical  p tr  i- 
meters m^crjb  mass,  mr*Lo.vl  mass , £*rope  length. 
Its  state  variables  are  xj*crab  position,  Xo*crab 
velocity,  x ^ "rope  angle  and  >:^*rope  angular 
velocity.  For  snail  rope  angles  the  linearized 
state  equations  are 
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with  r^*(rnc*HnL)g/nv;£ . Let  g^lQm/soc^  for  opera- 
tion on  earth.  Input  u is  the  force  accelerating 
the  crab.  Eigenvalues  ace  [ 0 ,0  , ji. , - j-aj]  . The 
observability  analysis  shows  that  x^  is  not  obser- 
vable by  *2 » x3  or  X/  * thus  the  crab  position 
sensor  is  essential. "*  It  was  shown  in  [ 5 1 that  \*o 
must  be  measured  or  estimated,  without  X2*feedback 
a stabilization  is  impossible. 

Given:  m^iOOOkg,  £**10  m,  maximum  load  3000  kg, 
design  a sampled-data  controller 

u(kl)  = -k'x(kT)  (43) 

for  the  following  specifications: 

i)  Consider  a typical  movement:  Pick  up  a load 
at  rest  and  drop  it  10m  away  at  rest  again, 
i.e.  initial  state  x*i 10  0 0 0]',  final  state 
x*0.  During  this  movement  the  required  force 
should  not  exceed  5000  Newton. 

ii)  The  amplitude  of  the  load  oscillation  after 
10  seconds  should  be  small  for  two  typical 
loads  of  3000  kg  and  1000  kg. 

iii)  It  is  desirable  to  avoid  the  measurement  of 
the  rope  angular  velocity  >:^ . 

The  sampling  interval  T was  selected  such  that  in 
the  worst  case  m ^=2 000  kg  the  complex  poles  in  z- 
plane  lie  on  a 45°  angle  with  respect  to  the  posi= 
tive  real  axis.  This  results  in  Tar,/8. 


The  discretization  and  evaluation  of  the  E matrix 
were  done  in  center  of  mass  coordinates,  in  which 
the  system  is  block  diagonalized.  The  result, 
transformed  back  to  sensor  coordinates  x,  is  then 
E-t^  <*E^+E.,]  a ■ m^£/ (mc-hnL)  (44) 


-3T/2 

-T/2 

T/2 

3T/2 

5T/2 


(nc*V  (-  sin  uT-  sin  2iT) 

“l  • (5  sin  x-T-4  sin  2xT+  sin  3u7) 

l 

2siai-t/2 (5si(K-T-Hsin2xT+sinii.r)‘ ' 
cos  (i-l/2)xI-2cos  (i-3/2)  iT+cos  (t-5/2  )xT)l 
-i  [_  ju[  sin(i-l/2)aT-2sia(i-3/2)  uX+sin(i-5/2)u)T)J 


1-0, 1,2. 3,4. 


Note  that  this  form  E=E(mr  ,mc ,T)  could  also  bo 
used  to  implement  a gain  scheduled  control  law 
k '"o 'E*k ' ,£ ,T) , which  keeps  the  eigenvalues 
constant . 


In  the  first  design  stop  a first  guess  for  k ' is 
determined  for  mL"J000kg  onLy.  First  a partial 
pol£  placement  is  made  by  eq.  (14),  which  gives  a 
!/•  2 damping  to  the  pendulum  motion  without 
changing  its  naturaL  frequency,  i.e.  i pole  pair 
at  zj  2=0. 4876^0. 3026.  For  the  initial  condition 
x(0)-f  10  0 0 0)'  the  first  control  input  is 

u(0)*l0k^,  thus  kj_  ^ 300  (N'ewton/m)  is  necessary  to 
meet  specif ication  i).  For  i fast  response  k^*500 
is  chosen.  After  specifying  two  eigenvalues  and 
one  feedback  gain  chore  remains  one  free  parameter, 
which  is  conveniently  exhibited  as  i parameter  on 
the  root  locus  for  the  remaining  two  poles.  Its 
complex  part  is  a circle  around  t*1  with  radius 
0.1786,  the  intersection  of  this  circle  with  the 
l/*.  2 damping  spiral  it  r.y  ■ *C. 3637^j0. 1177  is 
chosen.  This  results  in  *'=1300  1927  786"  -738). 


the  simulation  shows  that  u(kT)  does  not  »xcecd 
3000  Newton.  For  the  nominal  load  of  3000  kg  ti»e 
rn.axim.il  amplitude  after  10  seconds  is  4 . 3!',  of  the 
initial  displacement.  This  first  solution  is  how- 
ever unsatisfactory  for  n.  = 1000  kg,  with  a maximal 
amplitude  of  12.37,  if  ter  10  seconds. 

In  the  second  design  step  primarily  the  solution 
for  1 000  kg  must  be  speeded  up.  From  the  first 

solution  only  the  values  k^»500  and  k0-1927  are 
kept  and  k-j  and  are  Che  free  Darameters  of  the 
second  step.  The  four  eigenvalues  move  with  k->  and 
k^ , they  shall  be  kept  however  in  the  circle  r*0.5 
in  Fig.  2.  The  circle  maps  into  the  3000  kg 
boundary  in  the  k3«k^  plane  shown  in  Fig.  6. 

«4 


Fig.  6.  A crane:  Robustness  with  respect  to  large 
load  variations. 

At  point  A there  are  two  different  complex  con- 
jugate eigenvalue  pairs  crossing  the  circle  simul- 
taneously, such  that  the  complex  root  boundary 
intersects  itself.  The  right  real  boundary  is  out- 
side Che  figure.  AC  point  B the  complex  and  the 
left  real  root  boundary  for  z,*-0.05  meet,  i.e. 
this  k3.lv,  pair  generates  a double  pole  at  z=-0.05. 
In  points  C and  D a complex  pair  and  a real  root  at 
z**-0.05  cross  the  boundary  simultaneously.  The 
1000  kg  boundary  has  a similar  shape,  only  its 
right  part  and  the  intersection  D',  corresponding 
to  D,  are  shown.  Thus  the  first  design  point  indi- 
cated by  the  triangle  must  be  moved  to  the  left; 
let  ^*4000.  Since  k^  is  small  anyway,  k^*0  is 
chosen  in  view  of  specification  iii).  Simulation 
shows  that  the  maximum  amplitudes  after  10  seconds 
are  6.47,  for  3000  kg  and  37.  for  1000  kg  . The  angle 
x-j  remains  small,  such  that  this  assumption  for  the 
linearization  of  the  plant  equation  is  satisfied. 

If  necessary  a third  design  step  could  fallow  in 
which  k^*500  and  k^*0  are  fixed  and  k n and  are 
varied . 

5.  CONCLUSIONS 

Classical  parameter  piano  ideas  have  been  combined 
with  pole  placement  results  to  a design  method  in 
■V  space.  The  crucial  step  is  the  introduction  of  a 
plant  representation  in  the  form  of  the  matrix  E in 
eqs.  (16)  and  (A. 5).  The  linear  mapping  from  •? 
space  to  V space  is  performed  by  a multiplication 
with  the  matrix  E and  is  thus  reduced  to  a computa- 
tionally very  simple  step.  In  V space  typical 
design  aspects  such  as  actuator  limitations  and 
robustness  with  respect  to  sensor  and  actuator 
failures,  large  parameter  variations  and  short 
word  length  imp  lamentation  have  elementary  geometric 
interpretations,  and  several  questions  of  practical 
interest  can  be  treated  in  a clear  and  simple  way 
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is  is  illustrated  by  throe  examples. 

The  examples  ire  restricted  so  far  to  tradeoffs  in 
two  free  parameters  \t  a time,  where  a graphical 
interpretation  in  cross-sections  of  the  '*  space  is 
possible.  This  is  already  a practically  applic- 
able tool  with  apparent  advantages  over  graphical 
one-parameter  methods  like  root  locus.  For  example 
in  successively  closing  loops  jf  3 cascaded  system 
it  allows  to  make  tradeoffs  between  two  successive 
steps . 

The  concept  of  the  method  is  however  not  limited 
in  the  number  of  parameters.  Due  to  the  computa- 
tional simplicity  of  the  mapping  it  seems  feasible 
to  develop  computer-aided  design  methods  with 
displays  visualicing  three-dimensional  surfaces 
and  regions  by  moving  point  of  view  or  moving 
cross-section.  If  the  computer  has  to  make  the 
tradeoffs  in  problems  with  r.it.v  parameters , diffi- 
culties arise,  if  there  does  not  exist  a point  in 
V space  satisfying  all  speci fications . In  this 
situation  the  concept  of  a moving  boundary  may  be 
useful,  which  was  used  by  Zakian  and  AL-Naib  1 6] 
in  the  numerical  treatment  of  inequalities.  In 
Figs.  1 ar.d  Z this  means  that  the  pirameter  - or  r 
is  varied  continuously  until  a solution  is  found. 
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APPENDIX 

The  generalization  of  eq.  (16),  k ’ =2.’ E , to  the 
multivariable  case  was  published  in  (51  and  is 
quoted  here  for  easier  reference : 

Given  a controllable  pair  (A, 3),  let  ^“(b^  . . . t^J 
and  , , 

riJc=ranki  B,A  3 . . A*  * 3,  A*b  ...  A*%  . ] iai,2...p 

r.-r.  """  " <*•« 

ok  pk-l 

A controllability  index  i»l,2..p  of  the  pair 
(A,B)  is  the  smallest  integer  k such  that  r^j..a 

ri-ik’  il  b L is  linearly  dependent  5n  the 

vectors  left  of  it  in  the  controllability  matrix 
and  can  be  expressed  as 

A 3.  ..a""1  Blil-{A"ib1...A'lbl.l]st.  (A.:) 
In  order  co  aake  = ’ unique,  lec 
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Note  that  by  Popov's  theorem  on  feedback  invariants 
(71.  is  invariant  under  \ trans  format  Ion 
(A,B)  - (T(A-Bk[)T'1T  3),  deCT/0.  Let 
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Introduce  nxp  characteristic  parameters  in  a pxp 
matrix 
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?'  generalizes  o',  the  vector  of  coefficients  of  the 
characteristic  polynomial,  its  coefficients 
are  the  coordinates  of  an  nxp  dimensional  para^ 
meter  space  ■? . ?'  is  related  to  the  characteristic 

polynomial  by 

?(X)  “detCVI-A+B  K')»det(  P'-diagQ^)]  , 
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and  to  the  state  feedback  matrix  by 


K'  -MxB£'E.  (A. 3) 

Thus  the  system  representation  (E,MA3)  nay  be 
considered  as  a napping  between  two  nxp  dimen- 
sional parameter  spaces  r?  and  Y.  Note  chat  the  n 
coefficients  of  any  row  of  p*  enter  linearly  into 
the  determinant  in  eq.  (A. 7),  thus  they  can  be 
expressed  by  the  coefficients  of  P(\)  and  by  the 
remaining  nx(p-l)  characteristic  parameters  in  ? 1 , 
which  parameterize  the  remaining  degrees  of  free- 
dom after  pole  placement.  If  the  coefficients  of 
the  first  row  of  P1  are  eliminated  in  this  way, 
eq.  (A. 3)  results  in  n-p  equations  relating  the 
feedback  gains  to  the  n(p-l)  free  parameters. 

Due  to  the  structure  of  : ,A3 » n(p-l)  of  these 
equations  are  linear,  only  the  equations  for  the 
elements  of  the  first  row  of  K are  nonlinear  in 
the  free  parameters. 

The  free  parameters  can  now  be  chosen  according  to 
additional  requirements,  o.g.  minimizing  Che 
maximal  feedback  gain  in  view  of  actuator  limita- 
tions or  to  make  certain  columns  of  X1  equal  to 
zero  in  order  to  save  sensors. 


